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ABSTRACT 

We present a study of the stellar populations in the central ^ 200 pc of a large and 
homogeneous sample comprising 79 nearby galaxies, most of which are type 2 Seyferts. 
The star-formation history of these nuclei is reconstructed by means of state-of-the 
art population synthesis modeling of their spectra in the 3500-5200 A interval. A 
QSO-like featureless continuum (FC) is added to the models to account for possible 
scattered light from a hidden AGN. 

We find that: (1) The star- formation history of Seyfert 2 nuclei is remarkably het- 
erogeneous: young starbursts, intermediate age, and old stellar populations all appear 
in significant and widely varying proportions. (2) A significant fraction of the nuclei 
show a strong FC component, but this FC is not always an indication of a hidden AGN: 
it can also betray the presence of a young, dusty starburst. (3) We detect weak broad 
H/3 emission in several Seyfert 2s after cleaning the observed spectrum by subtracting 
the synthesis model. These are most likely the weak scattered lines from the hidden 
Broad Line Region envisaged in the unified model, given that in most of these cases 
independent spectropolarimetry data finds a hidden Seyfert 1. (4) The FC strengths 
obtained by the spectral decomposition are substantially larger for the Seyfert 2s which 
present evidence of broad lines, implying that the scattered non-stellar continuum is 
also detected. (5) There is no correlation between the star-formation in the nucleus 
and either the central or overall morphology of the parent galaxies. 

Key words: galaxies: active - galaxies: Seyfert - galaxies: stellar content - galaxies: 
statistics 



I INTRODUCTION 

Our understanding of how galaxies form and evolve has 
changed in a rather unexpected way in recent years. On 
the one hand, massive black-holes now appear to populate 
the nuclei of virtually every (massive enough) galaxy. On 
the other, the nuclei of active galaxies (AGN), previously 
thought to be the only galaxies that harbored massive black- 
holes, are now known to be the hosts of massive star-forming 
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regions. Hence, black-holes and starburst clusters coexist 
and are ubiquitous in the nuclear regions of galaxies. 

The observational evidence for this scenario is abun- 
dant. The presence of strong Call triplet absorptions in a 
large sample of Seyfert 2s and some Seyfert Is provided the 
first direct evidence for a population of red super-giant stars 
in their nuclear regions (Terlevich, Di'az & Terlevich 1990). 
The absence of signs of a Broad Line Region (BLR) in the 
direct optical spectra of Seyfert 2 nuclei which show broad 
lines in polarized light can only be understood if there is an 
additional central source of blue/UV continuum associated 
with the obscuring torus (Cid Fernandes & Terlevich 1995). 
This conclusion is further supported by observations that 
polarization is lower in the continuum than in the scattered 
broad lines (Miller & Goodrich, 1990; Tran 1995a,b) 

Most of the UV to near-IR continuum Ught in a num- 
ber of Seyfert 2s, in particular all those so far observed with 
HST, is due to nuclear starbursts that are resolved in HST 
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images. In some of these galaxies, the starburst luminosi- 
ties are comparable to those of their ionizing engines. Thus, 
an observer situated along the axis of the torus will detect 
comparable contributions to the optical continuum coming 
from the starburst and from the broad line region compo- 
nents (Heckman etal. 1995, 1997). Massive star forming re- 
gions are even found in the nuclei of powerful radio galaxies 
through the detection of Balmer absorption lines or strong 
Call triplet lines (Melnick etal. 1997; Aretxaga etal. 2001; 
Wills etal. 2002), showing that the intimate link between 
nuclear activity and star-formation exists not only in spiral 
galaxies, but also in massive ellipticals. 

We have been engaged in a large spectroscopic sur- 
vey of the nearest Seyfert galaxies in the southern hemi- 
sphere aimed at characterizing the properties of nuclear 
star-forming regions as a first step to understand the link 
between nuclear starbursts and massive black holes. The 
first results of this project were reported in Joguet etal. 
(2001, hereafter JOl), where high S/N observations of about 
80 Seyfert 2 galaxies were presented. This first inspection 
of the data revealed that about 50% of Seyfert 2s exihibit 
clear spectral signatures of star-formation in the recent past 
( < 100 Myr), most notably high-order Balmer lines in ab- 
sorption. The other half of the sample is equally split be- 
tween nuclei with strong emission lines but weak absorption 
lines and nuclei with a rich absorption line spectrum and 
reddish continuum, characteristic of old stellar populations. 
These results are in line with those obtained with smaller 
samples (Schmitt, Storchi-Bergmann & Cid Fernandes 1999; 
Gonzalez-Delgado, Heckman & Leitherer 2001; Cid Fernan- 
des etal. 2001a, hereafter CFOl), and highlight the high fre- 
quency of star-formation in Seyfert 2s. 

Recently we have also started to detect high order 
Balmer lines in Seyfert Is which may indicate the presence of 
extremely compact and luminous starbursts. The fraction of 
Seyfert 1 nuclei with Balmer absorption is understandingly 
lower than in Seyfert 2s given the additional difficulty of 
disentangling the stellar absorption features from the broad 
emission lines in Seyfert Is (Torres-Papaqui et al. 2004, in 
preparation). 

This paper is the first of a series devoted to a quan- 
titative analysis of the stellar populations in the nuclei of 
nearby Seyfert galaxies. We present detailed population syn- 
thesis models for a sample of 79 galaxies studied by JOl. 
For each galaxy we derive the nuclear star-formation his- 
tory (SFH), and correlate the nuclear stellar populations 
with other properties such as the emission line spectra and 
the (central) morphology of the parent galaxy. The models 
are also used to produce continuum-subtracted "pure emis- 
sion line" spectra which will be used in subsequent papers to 
study other properties of the nuclei including emission-line 
profiles, the relation between stellar and nebular kinemat- 
ics, etc. The ultimate aim of this study is to have at least a 
first glimpse at the coveted evolutionary connection between 
monsters and starbursts. 

This first paper concentrates on the presentation of our 
synthesis method and the results of its application to the JOl 
sample. Section|5|reviews the data set used in this study. A 
spectral synthesis method which combines tools from empir- 
ical population synthesis with the new Bruzual & Chariot 
(2003, hereafter BC03) models is described in Sectionl^and 
its application to 65 Seyfert 2s and 14 other galaxies is pre- 



sented in Section |1] The model subtracted spectra are used 
to search for weak broad emission features in Section|K| Sec- 
tion |H] presents an empirical characterization of the stellar 
populations in the sample based on a set of spectral indices. 
The statistics of star-formation in Seyfert 2s and relations to 
host morphology are discussed in Section |7| Finally, Section 
|S| summarizes our conclusions. 



2 THE DATA SET 

Highly efficient blue spectroscopic observations of nearby 
Seyfert 2 galaxies were carried out in four observing runs 
between May 1998 and August 2000, with the 1.5m ESQ 
telescope at La Silla (Chile) . The spectra cover a wavelength 
range from 3470 to 5450 A, with a spatial sampling of 0.82" 
per pixel. A detailed description of the acquisition and re- 
duction of the data, as well as of the selection of the sample 
has been published in JOl. 

For this work, ID spectra have been extracted afresh re- 
stricting ourselves to the central pixel, which, at the redshifts 
of the galaxy sample and adopting Ho = 75 kms~^Mpc~^, 
represents the central 59 to 439 pc for the Seyfert 2 galax- 
ies (median — 174 pc), and 45 to 596 pc for the Seyfert 
1, LINER, starburst and normal ones (median = 121 pc). 
Table [U lists some properties of the galaxies in the sample, 
including morphology, far-infrared properties and the linear 
scale. 

JOl listed 15 galaxies of the sample with dubious (or 
wrong) Seyfert 2 classification. After inspection of the spec- 
tra and a revision of literature data, six extra objects were 
reclassified. NGC 1097 (Storchi-Bergmann etal. 2003), NGC 
4303 (Colina etal. 2003; Jimenez-Bailon etal. 2003) and 
NGC 4602 were classified as LINERs. NGC 2935 and the 
mergers NGC 1487 and NGC 3256 were reclassified as Star- 
bursts (Dessauges-Zavadsky et al. 2000; Lira et al. 2002). Af- 
ter these minor revisions, the final sample has 65 Seyfert 2s, 
1 Seyfert 1, 4 LINERs, 4 Starburst/HII nuclei and 5 nor- 
mal galaxies. Note that our list of Seyfert 2s includes some 
mixed cases, such as NGC 6221 (Levenson et al. 2001), NGC 
7496 (JOl) and 7679 (Gu etal. 2001; Delia Ceca etal. 2001). 
Except for subtle details like high-excitation components in 
the wings of the [OIII] and H/3 profiles, these nuclei look 
like starbursts in the optical but have X-ray properties of 
AGN (see also Maiolino etal. 2003; Gongalves, Veron-Cetty 
& Veron 1999). 



3 SPECTRAL SYNTHESIS 
3.1 The method 

The modeling of stellar populations in galaxies has recently 
undergone a major improvement with the publication of 
high spectral resolution evolutionary synthesis models by 
BC03. We have incorporated this new library into a mod- 
ified version of the population synthesis code described by 
Cid Fernandes etal. (2001b). The code searches for the lin- 
ear combination of A'^* Simple Stellar Populations (SSP), ie, 
populations of same age (t) and metallicity (Z), which best 
matches a given observed spectrum O^. The equation for a 
model spectrum M\ is: 
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Sample Properties 



Name 


Type 


Hubble 


I-type 


/II 
pc/" 


a;(25,60) 


a(60,100) 


log Lpjff^ 


ESO 103-G35 


Sy2 


SO 


E 


257 


0.02 


1.54 


9.87 


ESO 104-Gll 


Sy2 


SBb 




293 


-2.00 


-1.22 


10.43 


ESO 137-G34 


Sy2 


SABO/a 


SO 


178 


-1.29 


-4.31 


10.12 


ESO 138-GOl 


Sy2 


E 


L 


177 


-0.47 


0.20 


9.66 


ESO 269-G12 


Sy2 


SO 




320 








ESO 323-G32 


Sy2 


SABO 




310 


-1.35 


-1.11 


9.81 


ESO 362-G08 


Sy2 


SO 


Sa 


309 


-1.45 


-1.67 


9.68 


ESO 373-G29 


Sy2 


SBab 


Sb 


181 








ESO 381-G08 


Sy2 


SBbc 


SBb 


212 


-1.30 


-0.85 


10.33 


ESO 383-G18 


Sy2 


SO/a"^ 




240 


-0.60 


-0.81 


9.40 


ESO 428-G14 


Sy2 


so/a 




105 


-1.04 


-0.62 


9.49 


ESO 434-G40 


Sy2 


SO/a 




160 








Fair 0334 


Sy2 


SB 


SBb 


359 


-1.15 


-1.46 


9.84 


Fair 0341 


Sy2 


so/a 


Sa? 


312 


-1.05 


-0.69 


9.75 


IC 1657 


Sy2 


SBbc 




232 


-2.66 


-1.94 


10.10 


IC 2560 


Sy2 


SBb 




189 


-1.41 


-1.24 


9.92 


IC 5063 


Sy2 


SO 




220 


-0.36 


0.48 


10.14 


IC 5135 


Sy2 


Sa 




313 


-2.34 


-0.86 


11.03 


IRAS 11215-2806 


Sy2 


SO** 


SO 


262 


-0.59 


-0.76 


9.40 


MCG +01-27-020 


Sy2 






227 


-0.51 


-2.07 


9.24 


MCG -03-34-064 


Sy2 


SB 




321 


-0.83 


0. 14 


10.53 


Mrk 0897 


Sy2 


Scd 




510 


-1.95 


-1.00 


10.69 


Mrk 1210 


Sy2 


Sa** 


Sa 


261 


0.11 


0.73 


9.83 


Mrk 1370 


Sy2 


Sa** 


Sa 


476 








NGC 0424 


Sy2 


SBO/a 


Sb 


226 


-0.04 


-0.02 


9.72 


NGC 0788 


Sy2 


so/a 


SO 


264 


0.00 


-0.29 


9.32 


NGC 1068 


Sy2 


Sb 




73 


-0.83 


-0.47 


10.77 


NGC 1125 


Sy2 


SBO/a 


Sb/c 


212 


-1.58 


-0.34 


9.95 


NGC 1667 


Sy2 


SABc 


Sc 


294 


-2.48 


-1.77 


10.62 


NGC 1672 


Sy2 


SBb 




87 


-2.40 


-1.47 


10.28 


NGC 2110 


Sy2 


SABO 


Sa 


151 


-1.82 


-0.62 


9.77 


NGC 2979 


Sy2 


Sa 




176 


-2.05 


-1.77 


9.56 


NGC 2992 


Sy2 


Sa 




149 


-2. 18 


-2.43 


10.26 


NGC 3035 


Sy2 


SBbc 




281 


-1.47 


-1.78 


9.62 


NGC 3081 


Sy2 


SABO/a 


SBO/a 


154 








NGC 3281 


Sy2 


Sab 




207 


-1.10 


-0.18 


10.24 


NGC 3362 


Sy2 


SABc 


Sb 


536 








NGC 3393 


Sy2 


SBa 


Sa 


242 


-1.25 


-1.06 


9.96 


NGC 3660 


Sy2 


SBbc 




238 


-2.44 


-1.74 


9.93 


NGC 4388 


Sy2 


Sb 




163 


-1.23 


-1.12 


10.28 


NGC 4507 


Sy2 


SABb 


SBa/b 


229 


-1.29 


-0.44 


10.14 


NGC 4903 


Sy2 


SBc 




319 


-1. 12 


-2.23 


9.85 


NGC 4939 


Sy2 


She 


Sa 


201 


-1.91 


-2.61 


9.92 


NGC 4941 


Sy2 


SABab 




72 


-1.09 


-2.17 


8.80 


NGC 4968 


Sy2 


SABO/a 


Sa 


191 


-0.93 


-0.41 


9.72 


NGC 5135 


Sy2 


SBab 


Sc 


266 


-2.23 


-1.03 


10.91 


NGC 5252 


Sy2 


SO 


SO 


447 








NGC 5427 


Sy2 


Sc 


Sc 


169 


-2.55 


-2.90 


10.24 


NGC 5506 


Sy2 


Sa 




120 


-0.96 


-0.11 


9.85 


NGC 5643 


Sy2 


SABc 




77 


-1.91 


-1.32 


9.93 


NGC 5674 


Sy2 


SABc 


SBc 


483 


-1.87 


-1.86 


10.44 


NGC 5728 


Sy2 


SABa 




180 


-2.54 


-1.15 


10.27 


NGC 5953 


Sy2 


Sa 


Sc 


127 


-2.46 


-1.26 


10.06 


NGC 6221 


Sy2 


SBc 


Sd 


96 


-2.20 


-1.65 


10.42 


NGC 6300 


Sy2 


SBb 


Sd 


72 


-2. 13 


-1.77 


9.78 


NGC 6890 


Sy2 


Sb 




156 


-2.03 


-1.46 


9.85 


NGC 7172 


Sy2 


Sa 




168 


-2.30 


-1.50 


10.09 


NGC 7212 


Sy2 


S? 


Irr? 


516 


-1.51 


-1.04 


10.72 


NGC 7314 


Sy2 


SABbc 


Sd 


92 


-2.13 


-2.60 


9.51 


NGC 7496 


Sy2 


SBb 




107 


-1.90 


-1.20 


9.83 


NGC 7582 


Sy2 


SBab 




102 


-2.32 


-0.77 


10.52 


NGC 7590 


Sy2 


Sbc 


Sd 


103 






9.85 


NGC 7679 


Sy2 


SBO 




332 


-2.13 


-0.56 


10.71 


NGC 7682 


Sy2 


SBab 


SBO 


332 








NGC 7743 


Sy2 


SBO 


SO 


110 


-1.57 


-2.28 


8.95 


Mrk 0883 


Syl 


? 




727 


-1.75 


-0.22 


10.50 


NGC 1097 


LINER 


E 




82 


-2.39 


-1.27 


10.34 


NGC 4303 


LINER 


SABbc 




101 


-3.23 


-1.97 


10.31 


NGC 4602 


LINER 


SABbc 




164 


-2.47 


-1.86 


10.03 


NGC 7410 


LINER 


SBa 


E/SO 


113 


-2.49 


-2.80 


8.99 


ESQ 108-G17 


STB 


10** 




142 


-2.20 


-0.90 


9.26 


NGC 1487 


STB 


P 




55 


-2.89 


-1.32 


8.85 


NGC 2935 


STB 


SABb 




147 


-2.55 


-2.18 


9.85 


NGC 3256 


STB 


SBm 




177 


-1.95 


-0.52 


11.23 


NGC 2811 


Normal 


SBa 




153 


-1.51 


-2.67 


9.04 


NGC 3223 


Normal 


Sb 




187 


-2.49 


-2.67 


10.14 


NGC 3358 


Normal 


SABO/a 




193 








NGC 3379 


Normal 


E 




59 








NGC 4365 


Normal 


E 




80 









Table 1. Column (2): Spectral class. Column (3): Hubble type, taken from RC3 except for objects marked with whose morphology 
was taken from NED. Column (4); Inner Hubble type, taken from Malkan, Gorjian & Tam (1998). Column (5): Linear scale. Columns 
(6) and (7): IRAS colours between 25, 60 and 100 //m. Column (8): Far infrared luminosity, in Lq. 

(ii) r\ = 10^" '*'^"*^ ^•'^^0^ is the reddening term, with A\ 
as the extinction at wavelength A. Extinction is modeled 
as due to an uniform dust screen, parametrized by the V- 
band extinction Av- The Cardelli, Clayton & Mathis (1989) 
extinction law with Ry =3.1 is adopted. 

(iii) X is the population vector, whose components Xj {j = 
1 . . . Ni,) represent the fractional contribution of each SSP 
in the base to the total synthetic flux at Aq. These flux 
fractions can be converted to a mass fractions vector fl using 
the model light-to-mass ratios at Ao. 



Mx {x, ,Av,v^,a^)= Mx^ 



Lj=i 



®G(w*,a*)(l) 



where 



(i) 6,,. s Lr^(t,, 
the j^^ SSP normalized at Ao 



)/Lf^^{tj,Zj) is the spectrum of 



The L^^^(tj, Zj) spectra are 



taken from BC03 models computed with the STELIB library 
(Le Borgne et al. 2003), Padova tracks, and Chabrier (2003) 
mass function. 
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(iv) Mxg, the synthetic flux at the normahzation wave- 
length, plays the role of a scaling parameter. 

(v) G(v*, a*) is the line-of-sight stellar velocity distribution, 
modeled as a Gaussian centered at velocity and broad- 
ened by (T*. 

The match between model and observed spectra is eval- 
uated by 

X\x, Mx,,Av,v^,a.) = Y^ [(^^ ~Mx)w>.f (2) 

A=l 

where the weight spectrum wx is defined as the inverse of the 
noise in 0\. Regions containing emission lines or spurious 
features (bad pixels or sky residuals) are masked out by 
assigning w\ — Q. 

We normalize Ox by its median value (Oat) in a window 
containing Ao instead of normalizing it by Oaqi which is 
obviously more sensitive to noise than the median flux. The 
synthetic spectrum is thus modeled in units of Ojv, and we 
expect to obtain M\g ~ 1 in these units. 

The search for the best parameters is carried out with 
a simulated annealing method, which consists of a series of 
Nm likelihood-guided Metropolis tours through the param- 
eter space (see Cid Fernandes et al. 2001b for an illustrated 
discussion of the Metropolis method applied to the popula- 
tion synthesis problem). In each iteration, the weights wx 
are increased by a factor thus narrowing the peaks in 
the {x, Mxq, Av ,Vi,,ai,) hyper-surface. The step-size (e) 
in each parameter is concomitantly decreased and the num- 
ber of steps {Ns) is increased to insure that each parameter 
can in principle random-walk through all its allowed range 
of values (eg, a;^ = — > 1 for the x components). To speed 
up computations, the kinematical parameters and cr* are 
kept fixed within the Metropolis loops, and re-fit at the end 
of each loop by a simpler minimization algorithm. 

The output population vector usually contains many 
Xj = components. The code offers the possibility of re- 
fitting the data excluding these components, which in prin- 
ciple allows a finer search for a best fit. In practice, however, 
this turns out to have little effect upon the results reported 
here. We have also implemented a conjugate-gradients rou- 
tine which runs after each Metropolis loop, in order to reflne 
the search for a minimum x^, but this feature too turns out 
to be largely irrelevant. The code "clips" deviant points, as- 
signing zero weight to pixels which are more than 3 sigma 
away from the rms Ox — Mx residual flux. This feature, 
which obviously requires an initial estimate of Mx, is useful 
to mask weak emission lines or defects in Ox not eliminated 
by the original wx ~ mask. 

This same general method, with variations on numerical 
schemes, spectral base, extinction laws and other details has 
been used in several recent investigations (eg, Kauffmann 
etal. 2003; Tremonti 2003; Mayya etal. 2004; BC03). 

3.2 The spectral base and tests 

A key ingredient in the synthesis of stellar populations is the 
spectral base bj^x- Ideally, the elements of the base should 
span the range of spectral properties observed in the sample 
galaxies and provide enough resolution in age and metallicity 
to address the desired scientific questions. 



We used as a starting base the 30 SSPs used by 
Tremonti (2003) in her study of SDSS galaxies. These SSPs 
cover 10 ages, f = 5 x 10^ 2.5 x 10^ 10*, 2.9 x 10^ 6.4x 10^ 
9 x 10^ 1.4 X 10^, 2.5 X 10^ 5 x 10^ and 1.1 x 10^° yr, and 3 
metallicities, Z = 0.2, 1 and 2.5Zq. Tremonti also included 
9 other spectra in her base, corresponding to continuous 
star formation models of different ages, e-folding times and 
metallicities. Not surprisingly, we found that these spectra 
are very well modeled by combinations of the SSPs in the 
base. Their inclusion in the base therefore introduces un- 
necessary redundancies, so we opted not to include these 
components. Smaller bases, with these same ages but only 
one or two metallicities, were also explored. 

Extensive simulations were performed to evaluate the 
code's ability to recover input parameters, to calibrate its 
technical parameters (number of steps, step-size, etc.) and 
to investigate the effects of noise in the data. To emulate 
the modeling of real galaxy spectra discussed in ^ the sim- 
ulations were restricted to the 3500-5200 A interval, with 
masks around the wavelengths of strong emission lines. 

Tests with fake galaxies generated out of the base show 
that in the absence of noise the code is able to recover the 
input parameters to a high degree of accuracy. This is illus- 
trated in figures^-c, where we plot input against output Xj 
fractions for 3 selected ages and all 3 metallicities. The rms 
difference between input and output Xj in these simulations 
is typically ~ 0.02. 

This good agreement goes away when noise is added to 
the spectra. As found in other population synthesis methods 
(eg., Mayya et al. 2004), noise has the effect of inducing sub- 
stantial rearrangements among x components corresponding 
to spectrally similar populations. The reason for this is that 
noise washes away the differences between bj^x spectra which 
are similar due to intrinsic degeneracies of stellar popula- 
tions (like the age-metallicity degeneracy; Worthey 1994). 
The net effect is that, for practical purposes, the base be- 
comes linearly dependent, such that a given component is 
well reproduced by linear combinations of others, as exten- 
sively discussed by Cid Fernandes etal. (2001b). For realis- 
tic S/N ratios these rearrangements have disastrous effects 
upon the individual Xj fractions, as illustrated in figures^l- 
f. The test galaxies in these plots are the same ones used in 
the noise-free simulations (panels a-c), but with Gaussian 
noise of amplitude S/N = 20 added to their spectra. Each 
galaxy was modeled 50 times for different realizations of the 
noise, and we plot the mean output Xj with an error bar 
to indicate the rms dispersion among these set of perturbed 
spectra. Clearly, the individual Xj fractions are not well re- 
covered at all. Hence, while in principle it is desirable to use 
a large, fine graded base, the resulting SFH expressed in x 
cannot be trusted to the same level of detail. 

A strategy to circumvent this problem is to group the 
Xj fractions corresponding to spectrally similar components, 
which yields much more reliable results at the expense of a 
coarser description of the stellar population mixture. This 
reduction can be done in a number of ways (Pelat 1997; 
Moultaka & Pelat 2000; Chan, MitcheU, & Cram 2003). Here 
we use our simulations to define combinations of x compo- 
nents which are relatively immune to degeneracy effects. We 
define a reduced population vector with three age ranges: 
t < 25 Myr; 100 Myr < t < lA Gyr and t > 2.5 Gyr. 
The resulting population vector is denoted by {xy,xi,xo), 
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where Y , I and O correspond to Young, Intermediate and 
Old respectively. Note that the metallicity information is 
binned-over in this description. Figures ^-i show the in- 
put against output {xy,xi,xo) components for the same 
S/N — 20 simulations described above. The deviations ob- 
served in the individual Xj fractions are negligible in this 
condensed description. Finally, we note that extinction does 
not suffer from the degeneracies which plague the popula- 
tion vector. For S/N — 20 the output Ay matches the input 
value to within 0.04 mag rms. 



4 APPLICATION TO SEYFERT 2S 

The synthesis method described above was applied to the 
spectra presented in ^^in order to estimate the star forma- 
tion history in the central regions of Seyfert 2s. 

4.1 Setup 

For the application described in this paper we have used the 
following set-up: Nm = 10 Metropolis loops, with weights 
increasing geometrically from /,„ = 0.1 to 10 times their 
nominal values while the step sizes decrease from e = 0.1 to 
0.005 for both the x components, M\g and Ay- The num- 
ber of steps in each Metropolis run is computed as a func- 
tion of the number of parameters being changed and by the 
step-size: Ns = {N^, + l)e~^. Since we are modeling galactic 
nuclei, it is reasonable to neglect the contribution of metal 
poor populations. We thus opted to use a base of A'^* = 20 
components with the same 10 ages described above but re- 
duced to Z = Zq and 2.5Zq. This setup gives a total of 
~ 3 X 10® steps in the search for a best model. 

4.2 Pre-processing 

Prior to the synthesis all spectra were corrected for Galactic 
extinction using the law of Cardelli et al. (1989) and the Ab 
values from Schlegel, Finkbeiner & Davis (1998) as listed 
m NEDV The spectra were also re-binned to 1 A sampling 
and shifted to the rest-frame. The fitted values of u* should 
therefore be close to zero. 

We chose to normalize the base at Ao = 4020 A, while 
the observed spectra are normalized to the median flux be- 
tween 4010 and 4060 A. This facilitates comparison with 
previous studies which used a similar normalization (eg., 
Gonzalez Delgado et al. 2004) . A flat error spectrum is as- 
sumed. We estimate the error in Ox from the rms flux in 
the relatively clean window between 4760 and 4810 A. The 
signal-to-noise ratio varies between S/N ^ 12 and 67, with 
a median of 27 in this interval, and about half these values 
in the 4010-4060 A range. 

Masks of 20-30 A around obvious emission lines were 
constructed for each object individually. In sources like NGC 
3256 and ESQ 323-G32, these w\ = masks remove rel- 
atively small portions of the original spectra. In several 

1 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 



cases, however, the emission line spectrum is so rich that the 
unmasked spectrum contains relatively few stellar features 
apart from the continuum shape (eg, Mrk 1210, NGC 4507). 
It some cases it was also necessary to mask the A < 3650 
A region because of strong nebular emission in the Balmer 
continuum (eg, NGC 7212). The strongest stellar features 
which are less affected by emission are Ca II K and the G- 
band. We have increased the weight of these two features 
by multiplying wx by 2 in windows of ~ 40 A around their 
central wavelengths, although this is of little consequence 
for the results reported in this paper. 

4.3 Featureless Continuum component 

A non-stellar component, represented by a F^/ oc v~^'^ 
power-law, was added to the stellar base to represent the 
contribution of an AGN featureless continuum, a traditional 
ingredient in the spectral modeling of Seyfert galaxies since 
Koski (1978). The fractional contribution of this component 
to the flux at Aq is denoted by xfc- 

In the framework of the unified model (eg, Antonucci 
1993), in Seyfert 2s this FC, if present, must be scattered ra- 
diation from the hidden Seyfert 1 nucleus. Nevertheless, one 
must bear in mind that young starbursts can easily appear 
disguised as an AGN-looking continuum in optical spectra, 
a common problem faced in spectral synthesis of Seyfert 2s 
(Cid Fernandes & Terlevich 1995; Storchi-Bergmann et al. 
2000). The interpretation of xpc thus demands care. This 
issue is further discussed in H5.3I 

4.4 Spectral Synthesis Results 

Figures mini illustrate some of the spectral fits. The exam- 
ples were chosen to represent the variety of spectra and star- 
formation histories found in the sample. The figures show 
the observed and synthetic spectra, as well as the Ox — Mx 
"pure emission" residual spectrum. Emission lines, partic- 
ularly weak ones, appear much more clearly in the resid- 
ual spectrum. The starlight subtraction also enhances the 
Balmer lines, particularly in nuclei with a significant inter- 
mediate age stellar population (eg, ESO362-G08). A detailed 
analysis of the emission lines is postponed to a future com- 
munication. 

The SFH plots (right panels in figures 1 ^101 show flux 
[xj) and mass (/ij) fractions for all 10 ages spanned by the 
base, plus the FC component. Note that by concentrating on 
the age distribution we are ignoring the metallicity informa- 
tion in X. Even binning over Z, a description of the SFH in 
terms of 10 age bins is too detailed given the efi'ects of noise, 
intrinsic degeneracies of the synthesis process, and limited 
spectral coverage of the data which is even further limited 
by the masks around emission lines. Rearrangements of the 
Xj strengths among adjacent age bins in these figures would 
produce fits of comparable quality, as found in independent 
population synthesis studies (eg, Mayya etal. 2004). 

As discussed in tl3.2l a coarser but more robust descrip- 
tion of the SFH requires further binning of the age distribu- 
tions in figures I2HH1 In Table |21 we condense the population 
vector X onto just four components: an FC (xpc) and three 
stellar components representing young {t < 2.5 x 10^ yr), in- 
termediate age (10* <t< 1.4x10" yr), and old (t > 2.5x10" 
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Figure 1. Tests of the spectral synthesis method. Panels a— c compare input against output Xj fractions for fake galaxies generated out of 
the base and with no noise. Each panel contains individual Xj components corresponding to 3 different metallicities but same age. Panels 
d-f are equivalent to panels a— c, but for simulations with noise of amplitude S/N = 20 added to each spectrum. The output fractions 
are now represented by the mean Xj, and the error bars correspond to the dispersion among 50 different realizations of the noise for 
each galaxy. The bottom panels (g— i) present the population vector condensed into just three age groups, xy, xj, and xq, binned-over 
all metallicities. 



yr) populations, denoted by xy, xj and xo respectively. We 
note that the 25 Myr component is zero or close to it in 
practically all galaxies, so that xy is virtually identical to 
the 5 Myr component. Mass fractions are described into just 
three components: /iy, ^/ and ^o, since we do not associate 
a mass to the FC. 

The values of Av and cr* produced by the synthesis 
are also listed in Table |5| The velocity dispersion was cor- 



rected by the instrumental resolutions of both the JOl spec- 
tra [ai-ast = 62 kms"^) and the STELIB library (~ 86 
kms~^). Figure [TTI compares our estimates of a* with val- 
ues compiled from the literature (mostly Nelson & Whittle 
1995). The agreement is good, with a mean and rms differ- 
ence of 3±21 kms~^. In the median, our values are 12 kms~^ 
larger than those in the literature, possibly due to filling of 
stellar absorptions by weak emission lines not masked out in 
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Synthesis Results 



Galaxy ^FC f^Y I ^* ^\ Notes 



ESO 


103-G35 


3 


4 





93 





.00 





.00 


100 


00 





.36 


114 





.9 


4 


,3 




BSO 


104-Gll 


3 


35 


32 


30 





.07 


8 


.25 


91 


68 


1 


.30 


130 


1 


.5 


12 


.7 




ESO 


137-G34 


15 








85 





.00 





.00 


100 


00 





.18 


133 


1 


.6 


4 


.1 




BSO 


138-GOl 


4 


49 


35 


12 





.32 


48 


74 


50 


95 





.00 


80 


2 


.6 


2 


.8 


WR? 


ESO 


269-G12 





1 


11 


88 





.00 





.52 


99 


48 





.00 


161 


1 


.4 


6 


.5 




BSO 


323-G32 








29 


71 





.00 


1 


.48 


98 


52 





.37 


131 





.5 


4 


.1 




ESO 


362^G08 








96 


4 





.00 


30 


38 


69 


62 





.36 


154 


1 


.4 


3 


.0 




BSO 


373-G29 


27 


22 


29 


21 





.06 


6 


.97 


92 


97 





.03 


92 


1 


.6 


3 


.3 




ESO 


3S1-G08 


24 


3 


67 


6 





.01 


30 


08 


69 


91 





.00 


100 





.7 


2 


.1 


BLR?, HBLR 


BSO 


383-GlS 


31 


24 





45 





.06 





.00 


99 


94 





.00 


92 


1 


.5 


5 


.0 


BLR 


ESO 


428G14 





20 


47 


33 





.03 


13 


23 


86 


75 





.53 


120 





.6 


3 


.7 


n-HBLR 


BSO 


434G40 





10 


11 


79 





.00 





.62 


99 


38 





.39 


145 





.3 


3 


.9 


HBLR 


Fai 0334 





22 


47 


31 





.03 


14 


72 


85 


25 





.31 


104 


1 


.1 


7, 


.8 




Fai 0341 


13 





1 


86 





.00 





.23 


99 


77 





.26 


122 


1 


.7 


8, 


.4 




IC 1657 





6 


33 


61 





.01 


4 


.98 


95 


01 





.51 


143 


1 


.1 


8, 


.4 




IC 2560 





33 





67 





.02 





.00 


99 


98 





.47 


144 


1 


.4 


4 


.1 




IC 5063 


6 


16 





78 





.01 





.00 


99 


99 





.57 


182 


1 


.2 


5 


.5 


HBLR 


IC 5135 


35 


34 


31 








.74 


99 


26 





00 





.00 


143 


2 


.6 


2 


.0 




IRAS 


11215-2806 


2 


8 


33 


56 





.02 


9 


.54 


90 


45 





.55 


98 


1 


.7 


5 


.9 




MCG 


+ 01-27-020 





20 


64 


15 





.06 


16 


55 


S3 


39 





.70 


94 


1 


.1 


5 


.3 




MCG 


-03-34-064 


5 


26 





69 





.02 





.00 


99 


98 





.38 


155 


1 


.5 


3 


.8 


BLR, HBLR, WR 


MRK 


0897 


18 


19 


60 


3 





.20 


60 


56 


39 


24 





.00 


133 


1 


.8 


2 


.1 




MRK 


1210 


3 


53 


5 


39 





.07 





.89 


99 


05 


1 


.03 


114 


4 


.1 


5 


.9 


HBLR, WR 


MRK 


1370 








73 


27 





.00 


13 


17 


86 


83 





.47 


86 


1 


.6 


4 


.7 




NGC 


0424 


36 


13 





52 





.02 





.00 


99 


98 





.62 


143 


1 


.8 


4 


.2 


BLR, HBLR, WR 


NGC 


0788 





13 





87 





.01 





.00 


99 


99 





.00 


163 


1 


.3 


6 


.0 


HBLR 


NGC 


1068 


23 


24 





54 





.04 





.00 


99 


96 





.00 


144 


1 


.1 


3 


.1 


BLR, HBLR 


NGC 




10 





64 


26 





.00 


48 


92 


51 


08 





.04 


105 


1 


.6 


5 


.8 




NGC 


1667 





8 


26 


66 





.01 


2 


.54 


97 


45 





.58 


149 


1 


.6 


8, 


.9 


n-HBLR 


NGC 


1672 


9 


12 


68 


11 





.03 


24 


55 


75 


42 





.24 


97 





.7 


2 


.7 




NGC 


2110 





26 


7 


67 





.01 





.05 


99 


93 





.38 


242 


2 


.8 


9 


.2 




NGC 


2979 





1 


78 


21 





.00 


15 


70 


84, 


30 





.70 


112 


1 


.4 


5 


.5 




NGC 


2992 





6 


35 


59 





.00 


2 


.60 


97 


40 


1 


.24 


172 





.9 


6 


.3 


HBLR 


NGC 


3035 


9 


5 


18 


69 





.00 





.18 


99 


81 





.08 


161 





.8 


4 


.6 


BLR 


NGC 


3081 


21 


5 





74 





.00 





.00 


100 


00 





.06 


134 





.5 


3 


.0 


HBLR 


NGC 


3281 


2 


16 





82 





.01 





.00 


99 


99 





.67 


160 





.8 


5 


.6 


n-HBLR 


NGC 


3362 


12 





30 


58 





.00 


10 


39 


89 


61 





.21 


104 


1 


.2 


5 


.7 


n-HBLR 


NGC 


3393 


14 


4 





83 





.00 





.00 


100 


00 





.06 


157 





.4 


2 


.7 




NGC 


3660 


45 





28 


26 





.00 


3 


.98 


96 


02 





.06 


95 


1 


.0 


2 


.9 


BLR, n-HBLR 


NGC 


4388 


16 


13 





70 





.02 





.00 


99 


98 





.73 


111 


1 


.5 


4 


.3 


HBLR 


NGC 


4507 


24 


19 





57 





.02 





.00 


99 


98 





.00 


144 





.6 


2 


.2 


BLR?, HBLR, WR? 


NGC 


4903 





6 


11 


83 





.00 





.50 


99 


50 





.55 


200 


1 


.0 


6 


.8 




NGC 


4939 





15 





85 





.01 





.00 


99 


99 





.18 


155 





.7 


3 


.7 




NGC 


4941 


2 


1 


11 


86 





.00 





.35 


99 


65 





.16 


98 





.1 


2 


D 


n-HBLR 


NGC 


4968 


1 


20 


58 


22 





.04 


15 


17 


84, 


79 


1 


.04 


121 


1 


.5 


7 


.9 




NGC 


5135 


12 


57 


24 


6 





.17 


6 


.27 


93 


56 





.38 


143 


2 


.8 


2 


.8 


n-HBLR 


NGC 


5252 





14 





86 





.01 





.00 


99 


99 





.57 


209 





.9 


3 


.8 


HBLR 


NGC 


5427 





16 


12 


72 





.01 


1 


.63 


98 


35 





.45 


100 


1 


.5 


6 


.9 




NGC 


5506 


63 







33 





.00 


4 


.75 


95 


25 





.76 


98 


1 


.9 


3 


D 


BLR, HBLR 


NGC 


5643 


6 





88 


6 





.00 


43 


15 


56 


85 





.30 


93 





.7 


2 


.1 


n-HBLR 


NGC 


5674 


2 





36 


63 





.00 





.71 


99 


29 





.64 


129 


1 


.1 


10 


.0 




NGC 


5728 





20 


51 


28 





.02 


6 


.99 


92 


99 





.25 


155 





.5 


2 


.0 


n-HBLR 


NGC 


5953 


18 


1 


74 


7 





.00 


57. 


94 


42 


06 





.67 


93 





.9 


2 


.7 




NGC 


6221 


32 


28 


32 


8 





.14 


14 


43 


85 


43 





.49 


111 


1 


.7 


3 


.3 




NGC 


6300 








75 


25 





.00 


16 


03 


S3 


97 





.92 


100 





.5 


4 


.2 


n-HBLR 


NGC 


6890 


13 


6 


9 


72 





.01 


2 


.27 


97 


72 





.45 


109 





.7 


3 


.9 


n-HBLR 


NGC 


7172 





11 


13 


76 





.01 





.95 


99 


05 





.80 


190 


1 


.7 


11 


.5 


n-HBLR 


NGC 


7212 


35 








65 





.00 





.00 


100 


00 





.36 


168 


1 


.5 


3 


.2 


BLR, HBLR 


NGC 


7314 


27 







68 





.00 





.00 


100 


00 





.83 


60 


2 


.0 


16 


.6 


HBLR 


NGC 


7496 


27 


23 


44 


6 





.07 


9 


.11 


90 


82 





.00 


101 


1 


.9 


2 


.3 


n-HBLR 


NGC 


7582 


31 





69 








.00 


100 


00 





00 


1 


.23 


132 


1 


.9 


2 


.1 


n-HBLR 


NGC 


7590 


1 





49 


50 





.00 


7 


.63 


92 


38 





.70 


99 





.4 


3 


.0 


n-HBLR 


NGC 


7679 


9 





91 








.00 


81 


81 


IS 


19 





.42 


96 


2 


.7 


2 


.6 




NGC 


7682 


6 


8 





86 





.00 





.00 


100 


00 





.00 


152 





.9 


4 


.2 


HBLR 


NGC 


7743 





1 


77 


21 





.00 


23 


10 


76 


89 





.48 


95 





.6 


3 


.4 




MRK 


0883 


48 


16 


31 


5 





.14 


74 


38 


25 


47 





.69 


202 





.7 


3 


.6 





NGC 1097 


19 


31 


37 


12 


0.05 


4.10 


95.85 


0.24 


150 


2.2 


3.9 


NGC 4303 


11 


4 


59 


26 


0.00 


12.69 


87.30 


0.10 


82 


0.3 


2.8 


NGC 4602 





25 


49 


25 


0.04 


8.10 


91.87 


0.36 


94 


1.1 


5.3 


NGC 7410 


12 


57 


24 


6 


0.17 


6.33 


93.50 


0.54 


140 


2.7 


2.8 


BSO 108-G17 


13 


72 


15 





4.70 


95.30 


0.00 


0.00 


108 


10.6 


2.1 


NGC 1487 





96 





4 


0.75 


0.00 


99.25 


0.00 


145 


6.0 


5.4 


NGC 2935 





12 


38 


50 


0.01 


1.93 


98.07 


0.37 


143 


0.5 


3.2 


NGC 3256 


24 


45 


30 





1.00 


47.94 


51.07 


0.45 


128 


4.7 


2.5 


NGC 2811 








13 


87 


0.00 


0.08 


99.92 


0.28 


246 


0.6 


4.2 


NGC 3223 








12 


88 


0.00 


0.15 


99.85 


0.18 


170 


0.7 


5.7 


NGC 3358 








18 


82 


0.00 


0.56 


99.44 


0.21 


194 


0.7 


5.0 


NGC 3379 








14 


86 


0.00 


0.27 


99.73 


0.00 


217 


0.7 


3.8 


NGC 4365 








17 


83 


0.00 


0.13 


99.87 


0.29 


257 


0.7 


4.6 



Table 2. Columns 2-5: Population vector in the (xpc , , ,^0) description, in percentage of the flux at 4020 A. Columns 6—8: 
percentage mass fractions associated to Y, I and O components. Column 9: Extinction (in V-band magnitudes). Column 10: Velocity 
dispersion, in kms~^. Column 11: P^r synthesized A. Column 12: Average absolute difference between model and observed spectra, 
in percentage. Column 13: BLR = weak broad H/3; HBLR = hidden broad line region detected through spectropolarimetry; n-HBLR = 
hidden broad line region not detected in spectropolarimetry; WR = Wolf-Rayet features. Question marks indicate dubious detections. 



the fits. On the whole, however, this is a minor effect given 
that the uncertainty in cr* is typically 20 1cm s~^ both for 
our and the literature values. 

The quality of the fits can be measured by Xai which 
is the of equation ^ divided by the effective number 
of wavelengths used in the synthesis (ie, discounting the 
masked points). In most cases we obtain Xa ~ li indica- 
tive of a good fit. However, the value of Xa depends on the 



assumed noise amplitude and spectrum, as well as on extra 
weights given to special windows, so standard x^ statisti- 
cal goodness-of-fit diagnostics do not apply. An alternative, 
albeit rather informal, measure of the quality of the fits is 
given by the mean absolute percentage deviation between 
0\ and Ma for unmasked points, which we denote by Aa. 
Qualitatively, one expects this ratio to be of order of the 
noise-to-signal ratio. This expectation was confirmed by the 
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simulations, which yield Aa = 10.9±4.8% for S/N = 10 and 
4.6 ±0.4% for S/N = 20. In the data fits we typically obtain 
values of 2-5% for this figure of merit (Table I^J, which is 
indeed of order of {S/N)~^. 

The random uncertainties in the model parameters were 
estimated by means of Monte Carlo simulations, adding 
Gaussian noise with amplitude equal to the rms in the 4760- 
4810 A range and repeating the fits 100 times for each 
galaxy. The resulting dispersions in the individual Xj com- 
ponents range from 2 to 8%, with an average (y{xj) of 4%. 
As expected, the binned proportions are better determined 
than the individual ones, with typical one sigma uncertain- 
ties of 2, 2, 4, and 4% for xfc, xy, xj and xo respectively. 
These values should be regarded as order of magnitude esti- 
mates of the errors in synthesis. The detailed mapping of the 
error domain in parameter space requires a thorough inves- 
tigation of the full covariance matrix, a complex calculation 
given the high-dimensionality of the problem. 



4.5 The diversity of stellar populations in Seyfert 
2s 

The spectral synthesis analysis shows that Seyfert 2 nuclei 
are surrounded by virtually every type of stellar population, 
as can be seen by the substantial variations in the derived 
SFHs (figures l^fTUl and Table 0. 

Some nuclei, like ESQ 103-G35 (Fig|21l, NGC 4903, and 
NGC 4941 are dominated by old stars, with populations 
older than 2 Gyr accounting for over 80% of the light at 4020 
A, while ESQ 362-G08 (Fig|HJ, Mrk 1370, NGC 2979 and 
others have very strong intermediate age populations, with 
xi > 70%. Young starbursts are also ubiquitous. In Mrk 
1210, NGC 5135 and NGC 7410 they account for more than 
half of the A4020 flux, and in several other Seyfert 2s their 
contribution exceeds 20%. Strong FC components are also 
detected in several cases (eg, NGC 5506 and Mrk 883). In 
general, at least three of these four components are present 
with significant strengths {x > 10%) in any one galaxy. 

This diversity is further illustrated in Fig ll2l where we 
present a face-on projection of the Xy/FC + + xq = 1 
plane, with xy/FC = xy+xfc- This heterogeneity contrasts 
with the results for our small sample of normal galaxies, all 
five of which are dominated by old stars {xo ~ 85%), with 
a minor contribution from intermediate age populations. 

4.6 Comparison with previous results 

Several galaxies in our sample had their stellar popula- 
tions studied before. In this section we compare our results 
with those reported by CFOl, who combines observations of 
Seyfert 2s from several other papers (Heckman et al. 1997; 
Storchi-Bergmann et al. 1998; Schmitt et al. 1999; Gonzalez 
Delgado etal. 1998; 2001). There are 13 galaxies in common 
with CFOl. Their data are similar in wavelength range and 
quality to those in this paper. Furthermore (not coinciden- 
tally, of course), they follow a description of stellar popu- 
lation plus FC components based on three broad groups, 
similar to those defined in tl3.2l and 14.41 which facilitates 
comparison with our results. The main differences between 
these two studies resides in the method to analyze the stel- 
lar population mixture. Whereas we make use of the full 
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Figure 11. Comparison of the stellar velocity dispersions esti- 
mated by our synthesis method and values compiled from the 
literature. The identity line is traced by the dashed line, while 
dotted lines indicate the ±1 sigma dispersion. 




Figure 12. Results of the spectral synthesis condensed on a 
(xyiFC 1 ^o) description and plotted in a face on projection of 
the Xy/FC ''r xj + xq = 1 plane. Squares denote normal galaxies 
and the stars mark starburst galaxies. AGN (mostly Seyfert 2s) 
are plotted with circles. Filled circles indicate objects in which 
broad lines have been identified either in direct or polarized light. 
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Figure 2. Results of the spectral fits for ESO 103-G35. The top left panel shows the logarithm of the observed and the synthetic spectra 
(shifted down for clarity). The Ox — Mx residual spectrum is shown in the bottom left panel. Spectral regions actually used in the 
synthesis are plotted with a thick line, while masked regions are plotted with a thin line. Crosses indicate points excluded from the 
synthesis by a 3 sigma clipping algorithm. All spectra are shown in units of the normalization flux Ojv. Panels on the right show the 
population vector binned in the 10 ages of the base. The top panel corresponds to the population vector in flux fraction (x), normalized 
to Aq = 4020 A, while the corresponding mass fractions vector fl is shown in the bottom. The power-law component xpc is plotted with 
an (arbitrary) age of 10^'^ yr and marked by a diamond. No mass is associated to this component (ie, fipc = in the bottom-right 
panel). 



spectrum, the population synthesis performed by CFOl is 
based in a handful of spectral indices (those listed in Table 
I^J. Also, the stellar population base used by CFOl is that 
of Schmidt et al. (1991), built upon integrated spectra of 
star clusters whose age and metallicity were determined by 
Bica & AUoin (1986), whereas here we use the BC03 mod- 



els. Notwithstanding the different methodologies, as well as 
differences in the spectra themselves (eg, due to the larger 
apertures in CFOl), one expects to find a reasonable agree- 
ment between these two studies. 

Figure 1131 summarizes the results of this comparison. 
Stars are used to plot the xy/FC fraction of CFOl (called 
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Figure 3. As figure^ but for NGC 5506 



just "xfc" in that paper) against our estimate, after re- 
scaling it to Ao = 4861 A to match the normalization in 
CFOl. The agreement is very good, with a mean offset of 
4% and rms of 6% between the two studies. Crosses and cir- 
cles are used for the xint and xold components of CFOl, 
which represent ~ 100 Myr and > 1 Gyr age bins. While 
qualitatively similar, these two components do not compare 
well with our definitions of xi and xo, which are systemat- 
ically larger and smaller than xint and xold respectively. 
This is essentially due to the different age-grouping schemes 
employed. Binning our 100 and 290 Myr base components 
into x'j and all t > 640 Myr components into x'o yields a 



good correspondence with the xjmt and xqld of CFOl, with 
residuals similar to those obtained for Xy/fc^ as shown in 
figure lOl 

The two discrepant crosses in the bottom of figure lT^ are 
IC 5135 (= NGC 7130) and NGC 5135. Again, the disagree- 
ment is only apparent. Our spectral synthesis find null con- 
tributions of the 100 and 290 Myr components, but for NGC 
5135 it fits 27% of the continuum at 4861 A as due to the 640 
Myr component, while for IC 5135 18 and 21% are associated 
with 640 and 900 Myr-old bursts respectively. In our discrete 
base these components are adjacent in age to the 290 Myr 
one used as the upper limit for x'j in the plot, and hence 
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Figure 4. As figurcEl but for Mrk 883. 



could well be regrouped within an intermediate age com- 
ponent by a slight modification of the age range associated 
with x'l, removing the discrepancy. In short, these ambigu- 
ities arise primarily because of the different age resolutions 
employed in the two studies. Finally, we note that our results 
for these two galaxies compare well with those of Gonzalez 
Delgado etal. (2001), who performed a detailed spectral fit 
of the Balmer absorption lines characteristic of intermediate 
age populations using their own models (Gonzalez Delgado, 
Leitherer & Heckman 1999). For both IC 5135 and NGC 
5135 they find ~ 50, 40 and 10% contributions (at 4800 
A) of young, intermediate age and old populations respec- 



tively, where the intermediate age component is modeled as 
a combination of 200, 500 and 1000 Myr bursts. Our con- 
densed population vector (at 4861 A) for these galaxies are 
(xY/FC^iy^o) ~ (60,40,0) and (56,27,17) respectively. 

We thus conclude that the population synthesis results 
reported here are consistent to within Aa; ~ 10% with pre- 
vious results, which is quite remarkable given differences in 
data and analysis methods. A clear advantage of our ap- 
proach with respect to methods like those in CFOl is that 
we model the full spectrum, which, besides providing more 
constraints, yields a more detailed understanding of the dif- 
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Figure 5. As figure but for MGC -03-034. 



ferent components which make up a Seyfert 2 spectrum. This 
advantage is explored next. 



5 ANALYSIS OF THE 

STARLIGHT-SUBTRACTED SPECTRA 

The excellent fits of the starlight in our spectra provide a 
unique opportunity to investigate the presence of weak emis- 
sion features in the Ox — Mx residual spectra. In this section 
we explore this opportunity and examine the consequences 
for the interpretation of the synthesis results. 



5.1 Wolf Rayet stars 

The presence of WR starts, which produce the so called WR 
bump at ~ 4680 A, provides a clear indication of starburst 
activity which is complementary to the population synthesis. 
We have used the starlight subtracted spectra to investigate 
the presence of the WR bump in our sample galaxies. 

WR stars were detected in the Starburst/WR galaxy 
ESO108-G17 and in the Seyfert 2s MCG-03-34, Mrk 1210 
and NGC 424. Hints of a WR bump are also present in ESO 
138-GOl and NGC 4507, but these cases require confirma- 
tion either from deeper observations or detection of the 5812 
A feature (Schaerer, Contini & Kunth 1999). These results 
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Figure 6. As figureEl but for IC 2560. 



are presented as notes in Tabled which summarizes the re- 
sults of the spectral synthesis models. Galaxies that show 
WR signatures present important young stellar components 
(average xy of 30%), thus validating the results of the syn- 
thesis. Given the severe blending with nebular lines of Hell, 
ArlV, NelV and Felll in Seyfert 2s, it is probably futile to 
attempt any detailed modeling of the WR population. 



5.2 Broad emission lines 

A close inspection of the starlight subtracted spectra reveals 
a weak broad H/9 in several of our Seyfert 2s. This is the case 



of NGC 3035, classified as a Seyfert 1 in NED but reclassified 
as a Seyfert 2 by JOl on the basis of the absence of broad 
lines. The top spectrum in figure [T^ is the same one used 
by JOl. When we subtract our starlight model, however, a 
weak but clear broad component emerges in H/3 (bottom 
spectrum in figure ITlk '). 

At least seven other objects exhibit a similar weak broad 
component under H/3: ESQ 383-G18, MCG-03-34-064, NGC 
424, NGC 1068, NGC 3660 (figure ITlb'). NGC 5506 and 
NGC 7212, and 2 other objects present marginal evidence 
for a broad H/3. In principle, these nuclei could be reclas- 
sified as Seyfert 1.5 or 1.8, following the notation of Oster- 
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Figure 7. As figurcEl but for NGC 5135. 



brock (1984). However, the detection of weak broad lines is 
not per se reason enough to classify a nucleus as of type 1. 
This violation of the historical definition of Seyfert types by 
Kachikian & Weedman (1974) is allowed in the context of 
the unified model, which predicts that a weak BLR should 
be detected in Seyfert 2s (at least those with a scattering 
region). We have thus opted not to reclassify these objects, 
with the understanding that these weak features are proba- 
bly scattered lines. 

The detection of a BLR provides a useful constraint for 
the interpretation of the spectral synthesis results. Broad 
emission lines in AGN always come with an associated non- 



stellar FC, presumably from the nuclear accretion disk. This 
is true regardless of whether the nucleus and BLR are seen 
directly, as in bona fide Seyfert Is, or via scattering, as in 
NGC 1068, the prototypical Seyfert 2 (Antonucci & Miller 
1985). Hence, when the BLR is detected we can be sure that 
an FC component is also present. 

Column 13 of Table|21hsts the spectropolarimetry infor- 
mation available for galaxies in our sample (from the com- 
pilation by Gu & Huang 2002 plus new data from Lumsden, 
Alexander & Hough 2004) . The existence of an AGN FC can 
also be established a priori in Seyfert 2s where spectropo- 
larimetry reveals a type 1 spectrum, ie, those with a "Hidden 
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Figure 8. As figure^ but for ESO 362-G08. 



BLR". This is the case of 17 objects in our sample. In 7 of 
these we were able to identify a weak broad H/3 in our direct 
spectra, while in the remaining cases the BLR detected with 
the aid of polarimetry is just too faint to be discerned in our 
data. In these cases, we expect a correspondingly weak FC 
component. Spectropolarimetry data are available for other 
fourteen galaxies in our list, but no hidden Seyfert 1 was 
detected, either because of the absence of an effective mir- 
ror or observational limitations. A further alternative is that 
these are "genuine Seyfert 2s" , ie, AGN with no BLR (Tran 
2001, but see Lumsden etal. 2004). 

Intriguingly, Tran (2001) reports no BLR detection 



from his spectropolarimetry data on NGC 3660, while figure 
114b leaves no doubt as to the existence of this component 
in direct light. Variability could be the cause for this dis- 
crepancy, in which case this galaxy should be classified as 
a type 1 Seyfert, because the BLR is seen directly. A type 
transition of this sort has happened in NGC 7582 (Aretxaga 
et al. 1999), which suddenly changed spectrum from that of 
a Seyfert 2 to one with broad emission lines sometime in 
mid 1998. 

Table |21 lists the galaxies for which we have evidence 
for the presence of weak broad emission lines, either from 
spectropolarimetry data in the literature (HBLR) or from 
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Figure 9. As figureEl but for Mrk 1370. 



the present study (BLR) . There is a clear tendency for these 
"broad hne Seyfert 2s" (BLS2s) to also show an important 
FC component in the synthetic spectra. We thus conclude 
that the synthesis is able to identify a bona fide non-stellar 
continuum when one is present. Further confirmation of this 
fact comes from the synthesis of the Seyfert 1 nucleus Mrk 
883, which yields xfc = 48%. As shown in Paper II (Gu 
et al. , in preparation), there is also a tendency for BLS2s to 
have larger [OIII] luminosities. This is consistent with them 
having more powerful AGN, which, for a fixed scattering 
efficiency, implies a brighter, easier to detect scattered FC, 
in agreement with the results reported above. 



5.3 FC versus Young Stars 

The synthesis also reveals a strong FC component in several 
objects for which there is no evidence of a BLR. The question 
then arises of whether xfc corresponds truly to an AGN 
continuum or to a dusty starburst. 

It is hard to distinguish a F„ oc power-law from 

the spectrum of a reddened young starburst over the wave- 
length range spanned by our data. The main differences be- 
tween the spectrum of a 5 Myr SSP seen through Ay ~ 2-3 
mag of dust and this power-law lie in the Balmer absorption 
lines and in the blue side of the Balmer jump. Since these 
features are masked in the synthesis, the xfc strength can 
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Figure 10. As figureEl but for the WR galaxy ESO 108-G17. 



be attributed to either a dusty-starburst, a true AGN FC, or 
to a combination of these. This confusion is hkely to occur in 
our synthesis models given that we impose a common Av for 
aU base components while young starbursts are intrinsically 
rather dusty (Selman etal. 1999). 

To illustrate this issue, we note that we obtain xfc = 
13% for the WR galaxy ESO 108-G17 (Fig[T01l. Given that 
there is no indication of an AGN in this irregular galaxy, it 
is very likely that this FC component is actually associated 
with a dusty starburst. The same applies to the starbursting 
merger NGC 3256, for which we derive xfc ~ 24%. Note, 



however, that in both cases the synthesis identifies a domi- 
nant young population, with xy = 72 and 45% respectively. 

Not surprisingly, the starburst-FC degeneracy is more 
pronounced in galaxies known to harbour an active nucleus. 
For NGG 6221, for instance, we find xfc = 32%, xy = 28% 
and Av = 0.5, while from the detailed imaging analysis 
and modeling of the full spectral energy distribution carried 
out by Levenson et al. (2001) we know that this nucleus is 
dominated by heavily extincted {Av ~ 3) young star clus- 
ters within < 1" of the AGN, with just a few percent of 
the optical light originating in the AGN responsible for its 
Seyfert-like X-ray properties. Most of the large xfc obtained 
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Figure 13. Comparison of the population vector obtained in this 
work with those in Cid Fernandes et al. (2001a) for 13 Seyfert 
2s in common. Different symbols indicate different age ranges. 
The identity line is traced by the dashed line, while dotted lines 
indicate the ±10% range. 

by the synthesis is thus associated with this reddened star- 
burst. Another emblematic example is IC 5135 (NGC 7130), 
whose UV spectrum shows wind lines of massive stars and 
a logarithmic slope typical of dusty starbursts, with little or 
no sign of a power-law (Gonzalez Delgado et al. 1998), while 
our synthesis yields xpc ~ 35%, xy — 34% and Av — 0. 
Similar comments can be made about other galaxies in the 
sample (eg, NGC 5135 and NGC 7582). 

These examples show that dusty starbursts, when 
present, are detected as an FC component by our synthe- 
sis method. The fraction xfc should thus be regarded as 
an upper limit to the contribution of a true FC, while the 
light fraction associated with very recent star-formation is 
bracketed between xy and Xy/pc- While there is no way to 
break this spectral degeneracy without broader spectral cov- 
erage, the information on the presence of weak broad lines, 
collected in Table |5| helps tackling this issue. As discussed 
in tl5.2l whenever a BLR component is detected either in 
direct or polarized light we can be sure that a true FC is 
present. One thus expects BLS2s to have larger xpc val- 
ues than Seyfert 2s for which we have no such independent 
evidence for the presence of an AGN continuum. 

This expectation is confirmed by the synthesis. The four 
largest values of xfc among Seyfert 2s, for instance, are all 
found in systems which show weak broad H/3. The link be- 
tween large xpc and the presence of broad lines is illustrated 
in figure lT^ which shows that the fraction of BLS2s increases 
systematically with increasing xpc- This effect is also clear 
in the statistics of xpc, which assumes a median (mean) 
value of 21% (19%) among BLS2s but just 2% (7%) among 
the rest of the Seyfert 2s. As predicted in il5.2l selecting only 
the galaxies where we see the BLR in our direct spectra in- 
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Figure 14. (a) Top: Total spectrum of NGC 3035. Bottom: 
Starlight subtracted spectrum of the same galaxy, showing a 
broad H/3 component not evident in the total spectrum. The 
dashed line marks the zero flux level for the bottom spectrum, 
(b) As the top panel, but for NGC 3660. 



creases the median (mean) xfc to 31% (29%), while the ten 
nuclei where a hidden BLR appears in polarized spectra but 
not in our data have weaker FCs: 5% (8%). 

Overall, these results confirm the prediction of Cid Fer- 
nandes & Terlevich (1995), who estimated that a broad com- 
ponent in H/3 should become discernible whenever the scat- 
tered FC contributes with > 20% of the optical continuum. 
Whenever xpc exceeds 20% but no BLR is seen, this 
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Figure 15. Histogram of the FC fraction in Seyfert 2s. Filled 
areas correspond to galaxies where a true FC is known to exist 
from the detection of a BLR in either direct or polarized light. 
The fraction of these latter galaxies with respect to the whole 
sample increases with increasing Xpcy shown by the dashed 
line and right-hand scale. 

component most likely originates in a dusty starburst rather 
than in a non-stellar source. 



6 STELLAR INDICES 

In order to facilitate the comparison of some relevant aspects 
of our population synthesis models and published work, we 
have found it convenient to express our results in terms of 
specific indices (lines and colours). In this section we present 
a simple characterization of the stellar populations in Seyfert 
2s on the basis of a set of indices commonly used in the 
literature. 

6.1 Direct measurements 

In order to provide an empirical characterization of stel- 
lar populations in our sample galaxies, we have measured 
a set of stellar indices directly from the observed spectra. 
To facilitate the comparison with previous studies, we have 
employed the index definitions used by Cid Fernandes et al. 
(2004), which are ultimately based on the studies by Bica & 
AUoin (1986a,b) and Bica (1988) of star cluster and galaxy 
spectra. 

We have measured the equivalent widths of the Ca II K 
(Wk), CN (Wcn) and G bands (Wo), as well as the colours 
C366O and C4510, defined as ratios between the continuum at 
3660 and 4510 A to the continuum at 4020 A, respectively. 
All these indices are measured with respect to a pseudo con- 
tinuum, which we define automatically following the recipe 
in Cid Fernandes et al. (2004) , after verifying that it works 



well for the present sample. Table |21 lists the resulting stel- 
lar indices. In this paper we give particular emphasis to the 
K line. Previous studies have shown that Wk, which essen- 
tially measures the contrast between old and young popula- 
tions, provides a useful mono-parametric description of stel- 
lar populations in galaxies. For instance, Seyfert 2s with un- 
ambiguous signatures of recent star-formation (such as UV 
wind lines, the WR bump, or high order Balmer absorption 
lines) aU have K lines diluted to Wk < 10 A (CFOl). 

Figure I16h shows the Wk versus C366O diagram. This 
plot has the advantage of involving only direct measure- 
ments, which, in the case of the JOl spectra, is not possible 
for more popular indices. The solid and dotted curves show 
the evolution of instantaneous burst BC03 models for solar 
(solid, green line) and 2.5 solar metallicity (dotted, red). The 
region occupied by pure stellar populations in this diagram 
is delimited by these tracks and a mixing line joining the 
youngest and oldest models. Points located below the BC03 
models in the bottom-right of the plot are explained as a 
combination of observational errors and intrinsic reddening, 
both of which affect more the C366O colour than Wk- 

The dashed line in figure 1161 represents a mixing line, 
where a canonical AGN Fi, oc power- law is added to 

an old stellar population in different proportions, from a 
pure power-law at {Wk, C3660) ~ (0, 1.05) to a pure 10^" yr 
population at ~ (20, 0.7). The fact that few galaxies line up 
along this mixing line shows that a simple "elliptical galaxy" 
plus power-law model, adopted in many early AGN studies, 
does not apply to the bulk of Seyfert 2s. On the contrary, 
the spread of points in this diagram implies a substantial 
diversity in the central stellar populations of Seyfert 2s, in 
agreement with our synthesis results (eg, figure IT^ as well 
as with previous studies (Cid Fernandes, Storchi-Bergmann 
& Schmitt 1998; JOl; Serote Roos etal. 1998; Boisson etal. 
2000). 



6.2 Indirect measurements 

Unfortunately, other useful tracers of the star-formation his- 
tory such as the Dn (4000) and H(5a indices (Balogh et al. 
1999; Worthey & Otaviani 1997), cannot be measured di- 
rectly for many of our galaxies because of severe contami- 
nation by emission lines. An indirect measurement of these 
indices can be performed over the model starlight spectra de- 
scribed in Naturally, the information contained in these 
"semi-empirical" indices is identical to the one contained in 
the model spectra from which they are measured. They are 
nevertheless useful for stellar population diagnostics based 
on observable indices, such as the Z)„(4000) versus H(5a di- 
agram amply explored by Kauffmann etal. (2002, 2003) in 
their study of SDSS galaxies. 

Figure 116b shows Wk and C366O measured from the 
deredenned synthetic spectra. The main difference with re- 
spect to the observed version of this diagram (figure IT^ ') is 
in the bottom-right points, which move upwards (to within 
the region spanned by the BC03 models) due to the cor- 
rection for intrinsic reddening. Objects in which a BLR is 
detected either in direct or polarized light (ie, BLS2s, repre- 
sented by filled circles) occupy a characteristic region in this 
plot, close to the power-law -|- old stars mixing line. This is 
consistent with the existence of an FC implied by the de- 
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Figure 16. C366O versus Wk diagram for the JOl sample. Panel 
(a) shows the results for C366O and Wk as measured directly from 
the observed spectra, whereas in (b) these two indices are mea- 
sured from synthetic spectra constructed to match the observed 
stellar spectrum. In both panels the solid (green) and dotted (red) 
curves represent the evolution of BC03 models of instantaneous 
bursts with Zq and 2.5Z0 metallicity respectively. The approx- 
imate location of populations with ages logt = 6, 7 , 8 , 9 and 
10 yr are indicated by the labels in panel (b). The dashed (blue) 
line between (Wfc, C'seeo) = (0, 1.05) and ~ (20, 0.7) represents a 
mixing line, where a oc u~^'^ power-law is added to a 10^" yr 
stellar population; diamonds mark power-law fractions in steps of 
25%. The average error bar is indicated in the top-right in (a). 
Symbols as in figure [T^ 



Table 3. Spectral indices, computed with the definitions of Cid 
Fernandes etal. (2004). Colours were corrected for Galactic ex- 
tinction. 



tection of a BLR. This tendency is also evident, albeit less 
clear, in the observed Wx-Cseeo diagram (figure ITCh ). 

Figure 1171 shows the synthetic D„ (4000)-H5a diagram 
for our sample. As in figure 1161 the scatter of points in 
this plot refiects the wide variety of stellar populations of 
Seyfert 2s, ranging from systems dominated by young stars, 
to post-starbursts and older systems, as well as mixtures 
of these populations. Again, objects where an FC is known 
to be present from the detection of a BLR are systemati- 
cally off-set towards the power-law + old stars mixing line. 
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Figure 17. D„{4000) versus H5a diagram. As in Figure IT6t ). 
both indices were measured from the synthetic spectra. The num- 
bers around the BC03 SSP evolutionary sequences mark the log- 
arithm of the age (in yr). Symbols and lines as in figure [TSl 

Note, however, that even these objects cannot be entirely 
explained in terms of this simple two-components model. 
An illustrative example in this respect is Mrk 1210, with 
D„(4000) = 1.3 and RSa = 1.6 A. This galaxy has both 
a hidden BLR detected through spectropolarimetry (Tran, 
Miller & Kay 1992) and a very young starburst, respon- 
sible for its WR bump (Storchi Bergmann, Cid Fernan- 
des & Schmitt 1998; Joguet 2001). Other examples include 
Mrk 477 (Heckman etal. 1997; Tran etal. 1992), Mrk 463E 
(Gonzalez-Delgado etal. 2001; Miller & Goodrich 1990) and 
NGC 424 (Moran etal. 2000; Joguet 2001). 



7 DISCUSSION 

7.1 Statistics of star- formation in Seyfert 2s 

Figure [TSl shows histograms of the Xy/FC, xj and xo com- 
ponents for the 65 Seyfert 2s in the sample, as well as the 
distributions of Wk and log t, the flux weighted mean stellar 
age (computed according to its definition in Cid Fernandes, 
Rodrigues Lacerda & Leao 2003). These histograms illus- 
trate the frequency of recent star-formation in Seyfert 2s. 

There are several possible definitions of "recent star- 
formation" , either in terms of observed quantities or model 
parameters. For instance, 31% of the objects have Wk < 10 
A, a rough dividing line between composite starburst + 
Seyfert 2 systems and those dominated by older populations 
according to CFOl. This frequency increases to 51% using a 
Wk < 12 a cut. In terms of xy/FC^ we find that nearly half 
(46%) of the sample has xy/FC > 20%, while for about one- 
third this fraction exceeds 30%. Old stars, although present 
in all galaxies, only dominate the light {xo > 70%) in 31% 



of the cases. In terms logt, we find that 25 out of our 65 
galaxies (38%) have \ogt < 8.5, ie, mean ages smaller than 
300 Myr. In summary, depending on the criterion adopted, 
between 1/3 and 1/2 of Seyfert 2s have experienced signif- 
icant star-formation in the past few hundred Myr. These 
numbers are similar to those derived by JOl by means of a 
visual characterization of stellar populations for the same 
sample and with those obtained by Storchi-Bergmann et al. 
(2001) for a sample ~ half as big as ours. 

Naturally, one would like to compare the incidence of re- 
cent star-formation in Seyfert 2s to that in other types of nu- 
clei. A particularly interesting comparison to be made is that 
between Seyfert 2s and LINERs, given that these two classes 
were once thought to belong to the same family (Ferland & 
Netzer 1983; Halpern & Steiner 1993). We have too few LIN- 
ERs in our sample to perform this comparison, but we may 
borrow from the results of Cid Fernandes et al. (2004) and 
Gonzalez Delgado etal. (2004), who have recently surveyed 
the steUar populations of LINERs and LINER/HII Transi- 
tion Objects using spectroscopic data practically identical to 
the one used in this paper and a similar method of analysis. 
Their results are very different from ours. For instance, while 
we find that over 2/3 of Seyfert 2s have Xy^FC > 10%, most 
LINERs have xy/Fc < 10%, and are often dominated by 
old stars. Still according to these authors, a large fraction of 
Transition Objects harbour significant numbers of 10^-10^ 
yr stars, similar to those found in a;/-dominated Seyfert 2s 
(eg, ESO 362-G08 and Mrk 1370, figures El and . Their 
Xy/FC strength, however, is only 6% on average, compared 
to 23% for our Seyfert 2s. 

Thus, a substantial fraction of nearby Seyfert 2s live in 
much younger stellar environments than LINERs or Tran- 
sition Objects, which suggests an evolutionary sequence. 
Given that the later objects are known to harbour less lu- 
minous AGN and have different gas excitation conditions 
than Seyferts, if confirmed, this would imply a substantial 
evolution of the AGN itself — in parallel with the evolution 
of stars in its surroundings. 

7.2 Aperture effects 

Spectra collected through apertures of fixed angular size 
sample distance-dependent linear scales, which may intro- 
duce systematic effects on the derived stellar population 
mixtures. In order to investigate whether such potential bi- 
ases affect our results in some systematic way we plot in fig- 
ure ll9l the equivalent width of Call K against the pc/arcsec 
scale. (Recall that our observations cover an area of ^ 1 
arcsec^.) No systematic trend is observed between this stel- 
lar population tracer and the linear dimension of the regions 
sampled by our data. The same conclusion holds for the x 
components, none of which correlates with distance either. 
The absence of systematic trends indicates that our analysis 
is not affected by aperture effects. 

7.3 Relation to host galaxy morphology 

A galaxy's morphology carries important information for 
the study of star formation history in galaxies (eg., Kenni- 
cutt 1998). From early to late types, galaxy colours become 
bluer and star formation rates higher. Storchi-Bergmann 
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Figure 18. Distributions of the Xy/pc^ ^^'^ fractions, Wk and the flux weighted mean stellar age inferred from the synthesis. 
Only Seyfert 2s were included in these histograms. 



et al. (2001) studied the relationship between nuclear stel- 
lar population and morphology in a sample of 35 Seyfert 2s, 
13 in common with our sample. They found that the frac- 
tion of galaxies with recent nuclear/circumnuclear starburst 
increases along the Hubble sequence and that there is a re- 
lation between the presence of a nuclear starburst and a late 
"inner Hubble type" from the HST classifications of Malkan, 
Gorjian & Tam (1998). 

In order to check this result in our larger sample, we 
have collected morphological information on the host galax- 
ies from RC3, and the inner Hubble types from Malkan et al. 



(1998). In figure we plot three of our stellar population 
tracers (Wk, xy/fCj a-nd mean age) against the RC3 and 
inner morphological types (later types have larger numbers). 
There is no correlation in our sample between (optical) mor- 
phology and nuclear star formation. Since the HST snapshot 
survey of Malkan etal. (1998) wEis taken using the F606W 
filter which is sensitive to dust extinction and Ha emission, 
we searched the HST/NICMOS archive data and found im- 
ages for 36 of our galaxies taken with the F160W filter. No 
clear trend between morphological type and nuclear stellar 
population is found in the near-IR images either. 
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Figure 19. Projected linear scale against the equivalent width 
of Call K. Symbols as in figure [T^ 



Finally, following Storchi-Bergmann et al. we also 
looked for companions to our Seyfert 2s according to an- 
gular separation, radial velocity, and apparent magnitude. 
We find that 40% of the galaxies with W^k < 10 A and 47% 
of those with Wk > 10 A show the presence of companions, 
compared respectively with 42% and 47% that do not have 
near neighbours. 

Hence, contrary to the findings of Storchi-Bergmann 
et al. (2001) we find no correlation between the nuclear stel- 
lar populations of Seyfert 2 galaxies and either morpho- 
logical types or environments. Since there is a substantial 
overlap between our samples, the discrepancy is most likely 
statistical. Given the importance of the result, however, it 
would be worthwhile to repeat the study in an even larger 
sample of galaxies. 



8 SUMMARY 

We have presented a study of the stellar population in the 
central ~ 50-500 pc of a large, well defined and homoge- 
neous sample of Seyfert 2s from the atlas of JOl. Our main 
results may be summarized as follows: 

(i) We have developed a spectral synthesis code which 
decomposes an observed spectrum as a sum of simple stel- 
lar populations, represented by state-of-the-art evolutionary 
synthesis models, plus an AGN continuum. Unlike previous 
studies, which use only a few absorption equivalent widths 
and colours (eg Bica 1988; Schmitt et al. 1999; Cid Fernan- 
des etal. 2001a,b), we now synthesize the whole observed 
spectrum. The method produces an estimate of the star- 
formation history in terms of flux (x) or mass (/I) fractions 
associated to each component in the spectral base, as well as 
extinction and velocity dispersions. As in other population 



synthesis studies, the accuracy of our method is limited by 
noise, spectral coverage, and intrinsic degeneracies of stellar 
populations. These limitations are not critical as long as one 
does not attempt an over-detailed description of the SFIf. 

(ii) The synthesis method was applied to 3500-5200 A 
spectra of 65 Seyfert 2s and 14 other galaxies from the JOl 
atlas. The synthesis produces excellent fits of the observed 
spectra, with typical flux residuals of 2-5%. The star for- 
mation history deduced from these flts is remarkably varied 
among Seyfert 2s. Young starbursts, intermediate age, and 
old populations all appear in significant and mixed amounts. 

(iii) The stellar velocity dispersions obtained from the 
synthesis are in good agreement with measurements in the 
literature for galaxies in common. This is important since 
this information gives us a handle on the black-hole mass 
(Ferrarese & Merrit 2000; Gebhart etal. 2000). 

(iv) The starlight-subtracted spectra were used to inves- 
tigate the presence of weak broad emission features which 
are hard to detect in the total spectrum. This analysis re- 
vealed the signatures of WR stars in 3 (maybe 5) Seyfert 2s, 
all of which have large young stellar populations as deduced 
by the synthesis. 

(v) The analysis of the "pure-emission" spectra further 
allowed the detection of a weak BLR-like component under 
H/3 in several Seyfert 2s. For most of these objects indepen- 
dent spectropolarimetry data reveals a type 1 spectrum. We 
are thus detecting in direct, non-polarized light the scattered 
BLR predicted by the unified model. 

(vi) The detection of a BLR either in direct or polarized 
light implies that an associated FC should also be present 
in our spectrum. Indeed, the FC strengths obtained by the 
spectral decomposition for the 17 Seyfert 2s in our sample 
which present such evidence are substantially larger than 
for nuclei with no indications of a BLR. The fraction of 
these "broad line Seyfert 2s" increases systematically with 
increasing xpc- 

(vii) Dusty young starbursts can also appear disguised 
as an AGN-looking continuum. We have identified several 
cases where this indeed happens. These Seyfert 2s tend to 
have large xfc o-nd xy fractions. The information on the 
existence of a BLR, and thus of a true FC component, 
helps breaking, at least partly, this starburst-FC degeneracy. 
Whenever the spectral synthesis identifies a strong ( > 20%) 
FC but no BLR is detected either in the residual or polar- 
ized spectra, the FC is most likely dominated by a dusty 
starburst rather than scattering of a hidden AGN. 

(viii) Stellar indices were measured both from the ob- 
served and synthetic spectra to provide a more empirical 
characterization of the stellar content of our galaxies and 
facilitate comparisons with independent studies. These in- 
dices, which include Wk, the 3660/4020 colour, L>„(4000) 
and H(5a, confirm the heterogeneity of nuclear stellar popu- 
lations in Seyfert 2s. 

(ix) Between 1/3 and 1/2 of nearby Seyfert 2s have expe- 
rienced significant star formation in the recent past. Thus, 
as a class, Seyfert 2 nuclei live in a much younger stellar 
environment than LINERs and normal galactic nuclei. 

(x) We find no significant correlation between the host 
morphology (as deduced both from ground based and HST 
optical and near-IR images) and stellar population in Seyfert 
2s. The presence of companions does not seem to correlate 
with the stellar populations either. 
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Figure 20. Distribution of morphological Hubble type for the host galaxies (top panels) and the inner region (bottom panels) against 
stellar population properties; Wx (left), Xy/pc (middle) and the mean age (right). Symbols as in figure [T^ 



We have shown that using modern population synthe- 
sis models it is possible to obtain meaningful information 
about the history of star formation in objects for which the 
spectral data are dominated by strong emission lines. The 
road to detailed quantitative investigations of the relation- 
ship between star-formation and AGN is now open. We will 
be traveling along it in our forthcoming papers. 
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